home *** CD-ROM | disk | FTP | other *** search
/ Amiga Format CD 46 / Amiga Format CD46 (1999-10-20)(Future Publishing)(GB)[!][issue 1999-12].iso / -in_the_mag- / reader_requests / scilab / demos / robust / mu.dem < prev    next >
Text File  |  1999-09-16  |  1KB  |  41 lines
  1. //test of musolve
  2. mode(-1)
  3. M1 = [ ..
  4.    5.2829       5.7683      -2.4004       1.2205      -6.4148    
  5.    9.7769e-01   2.9786      -3.0408       5.0257e-01  -2.6504    
  6.    7.0819       9.6324      -3.5750       3.3016      -6.7030    
  7.   -1.6261      -2.9763       1.6870      -1.0603       1.2211    
  8.    2.3056       4.3712      -2.4785       2.6152      -1.9832    ];
  9.  
  10. M2 = [ ..
  11.   -1.1308      -1.7785       8.7974e-01  -7.5206e-01   1.2089    
  12.   -3.5255e-01  -5.7002e-01   2.9305e-01  -2.5442e-01   3.7691e-01
  13.   -1.3724      -2.1501       1.0741      -9.1188e-01   1.4669   
  14.    3.5839e-01   5.5101e-01  -2.7290e-01   2.3565e-01  -3.7663e-01
  15.   -4.9015e-01  -7.8706e-01   4.0215e-01  -3.3617e-01   5.3261e-01];
  16. //*******************************************************
  17. M=M1 +%i*M2;
  18. // Let the structure be all scalar blocks
  19. K = [1 1 1 1 1]'; 
  20. // Let the first, the third and the fifth blocks be real,
  21. // and let the rest of blocks be complex
  22. T = [1 2 1 2 1]';
  23. [D,g,mu] = musolve(M,K,T);
  24. spec(M'*D*M+%i*(G*M-M'*G')-mu^2*D)
  25.  
  26.  
  27. // Now, we compute it again with respect to all complex blocks
  28. T = [2 2 2 2 2]';
  29. [D,g,mu] = musolve(M,K,T);
  30. spec(M'*D*M+%i*(G*M-M'*G')-mu^2*D)
  31.  
  32.  
  33. T = 3*[1 1 1 1 1]';
  34. [D,g,mu] = musolve(M,K,T;
  35. //
  36. K = [2 3]';
  37. T = [2 2]';
  38. [D,g,mu] = musolve(M,K,T);
  39.  
  40.  
  41.